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ABSTRACT 

We study generalized log-sine integrals at special values. At 
TT and multiples thereof explicit evaluations are obtained in 
terms of Nielsen poly logarithms at ±1. For general argu- 
ments we present algorithmic evaluations involving Nielsen 
poly logarithms at related arguments. In particular, we con- 
sider log-sine integrals at 7r/3 which evaluate in terms of 
poly logarithms at the sixth root of unity. An implementa- 
tion of our results for the computer algebra systems Mathe- 
matica and SAGE is provided. 

Categories and Subject Descriptors 

1. 1.1 [Symbolic and Algebraic Manipulation]: Expres- 
sions and Their Representation; 1.1.2 [Symbolic and Al- 
gebraic Manipulation]: Algorithms 

General Terms 

Algorithms, Theory 



Keywords 

log-sine integrals, multiple poly logarithms, multiple zeta val- 
ues, Clausen functions 



1. INTRODUCTION 

For n = 1, 2, . . . and fc > 0, we consider the (generalized) 
log-sine integrals defined by 



Lsi"' (a) 



nk 1 n—\—k 

e log 



2sin- 



de. 



(1) 



The modulus is not needed for < cr < 27r. For k = 
these are the (basic) log-sine integrals Ls„ (a) := Ls$f' (a). 
Various log-sine integral evaluations may be found in [20, 
§7.6 & §7.9]. 

In this paper, we will be concerned with evaluations of 
the log-sine integrals Lsi*^' (a) for special values of a. Such 
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evaluations are useful for physics [15]: log-sine integrals ap- 
peared for instance in recent work on the e-expansion of 
various Feynman diagrams in the calculation of higher terms 
in the e-expansion, [8, 16, 9, 11, 14]. Of particular impor- 
tance are the log-sine integrals at the special values tt/3, 
tt/2, 27r/3, TT. The log-sine integrals also appear in many 
settings in number theory and analysis: classes of inverse 
binomial sums can be expressed in terms of generalized log- 
sine integrals, [10, 4]. 

In Section 3 we focus on evaluations of log-sine and related 
integrals at ir. General arguments are considered in Section 
5 with a focus on the case 7r/3 in Section 5.1. Imaginary 
arguments are briefly discussed in 5.2. The results obtained 
are suitable for implementation in a computer algebra sys- 
tem. Such an implementation is provided for Mathematica 
and SAGE, and is described in Section 7. This complements 
existing packages such as Isjk [15] for numerical evaluations 
of log-sine integrals or HPL [21] as well as [25] for working 
with multiple polylogarithms. 

Further motivation for such evaluations was sparked by 
our recent study [6] of certain multiple Mahler measures. For 
k functions (typically Laurent polynomials) in n variables 
the multiple Mahler measure /i(fi, P2, ■ ■ ■ , Pk), introduced 
in [18], is defined by 



Jo io ' ^ ' ' 



. dtn- 



When P = Pi = P2 ~ ■ ■ ■ = Pk this devolves to a higher 
Mahler measure, fik{P), as introduced and examined in [18]. 
When k — 1 both reduce to the standard (logarithmic) 
Mahler measure [7]. 

The multiple Mahler measure 

/ifc(l + a: + y*) ■— fj,{l + x + yi,l + x + y2, . . . ,l + x + yk) (2) 

was studied by Sasaki [24, Lemma 1] who provided an eval- 
uation of jj,2{l + X + y,). It was observed in [6] that 

/ife(l + x + y,) ^ - Lsfc+i (^) - - Lsfc+i (tf) . (3) 

Many other Mahler measures studied in [6, 1] were shown 
to have evaluations involving generalized log-sine integrals 
at TT and 7r/3 as well. 

To our knowledge, this is the most exacting such study 
undertaken — perhaps because it would be quite impossible 
without modern computational tools and absent a use of the 
quite recent understanding of multiple polylogarithms and 
multiple zeta values [3]. 



2. PRELIMINARIES 

In the following, we will denote the multiple polylogarithm 
as studied for instance in [4] and [2, Ch. 3] by 



For our purposes, the ai , . . . , will usually be positive in- 
tegers and ai > 2 so that the sum converges for all jzj < 1. 

For example, Li2,i(2) = X^feLi 'p' Ej=i j- particular, 
Life(a;) :— X^S^i ^ the polylogarithm of order k. The 
usual notation will be used for repetitions so that, for in- 
stance, Li2_{i}3(2) = Li2,i,i,i(2:)- 

Moreover, multiple zeta values axe denoted by 



C(ai, 



,afc) — Liai,...,afc(l). 



Similarly, we consider the multiple Clausen functions (CI) 
and multiple Glaisher functions (Gl) of depth k and weight 
w = a\ + . . . + Ok defined as 



Clai,...,afc {0) — 
Glai,...,afc {0) = 



Im Liai,. 




if 


w 


even 


Re Liai,. 




if 


w 


odd 


Re Liaj,. 




if 


w 


even 


Im Liai,. 




if 


w 


odd 



(4) 
(5) 



in accordance with [20]. Of particular importance will be 
the case of 9 — n /3 which has also been considered in [4]. 

Our other notation and usage is largely consistent with 
that in [20] and that in the newly published [23] in which 
most of the requisite material is described. Finally, a recent 
elaboration of what is meant when we speak about evalua- 
tions and "closed forms" is to be found in [5]. 

3. EVALUATIONS AT tt 

3.1 Basic log-sine integrals at n 

The exponential generating function, [19, 20], 



^ oo 

— 7^ LSm+l (tt) 



r(i + A) 
r2(i + A) 



(6) 



(7) 



is well-known and implies the recurrence 
(-1)" 

; — Ls„+2 (vr) = TT a(n + 1) 

n! 

1-2 {_\\^ 

+ 1] (fcTT)! ""^^ ~ ^""^ ' 

where a(m) = (1 - 2^-'")C(m). 

Example 1. (Values of Ls„ (tt)) We have Ls2 (vr) = and 

- LS3 



LS4 
-LSg 

LS6 
- LS7 

Lsg 



12 

= |^C(3) 



19 5 

= 240 ^ 

= Y^C(5) + |7r\(3) 

275 7^45 2 
= 1344 " +T"<(3) 

-^-C(7) + ^^\(5)-f^.^(3), 



and so forth. The fact that each integral is a multivari- 
able rational polynomial in tt and zeta values follows directly 
from the recursion (7). Alternatively, these values may be 
conveniently obtained from (6) by a computer algebra sys- 
tem. For instance, in Mathematica the code 
FullSimplif y[D [-Binomial [x, x/2] , {x.e}] /.x->0] 
produces the above evaluation of Lse (tt). O 

3.2 The log-sine-cosine integrals 

The log-sine-cosine integrals 

LsCm.n (o") ^ log*"^ 2 siu — log"^'^ 2 cos— dO 
Jo 2 2 

(8) 

appear in physical applications as well, see for instance [9, 
14]. They have also been considered by Lewin, [19, 20], and 
he demonstrates how their values at cr = tt may be obtained 
much the same as those of the log-sine integrals in Section 
3.1. As observed in [1], Lewin's result can be put in the form 



1 

TT 



Lsc 



m-t- 1 ,n4- 1 



2"+T r (i±s) 



r 1 



x+y ) 
2 ) 



X \l y 

x/2 U/2 



r(i + f)r(i + f) 
r(i + 4^) 



(9) 



The last form makes it clear that this is an extension of (6). 

The notation Lsc has been introduced in [9] where evalu- 
ations for other values of a and low weight can be found. 

3.3 Log-sine integrals at tt 

As Lewin [20, §7.9] sketches, at least for small values of 
n and fc, the generalized log-sine integrals Lsi*"' (tt) have 
closed forms involving zeta values and Kummer-type con- 
stants such as Li4(l/2). This will be made more precise in 
Remark 1. Our analysis starts with the generating function 
identity 



E Ls^.+i (-) 



n,k>0 

- ie 



n\ k\ 



2 sin- 



i(.LO 



AO 



7 Bi (/i - |, 1 + A) - ie'""Bi/2 (a* - i - 



(10) 



given in [20]. Here Bx is the incomplete Beta function: 



Bx{a,b) 



dt. 



We shall show that with care — because of the singularities 
at zero — (10) can be differentiated as needed as suggested 
by Lewin. 

Using the identities, valid for a,b > and < a:: < 1, 



B4a,b) = LAl-^ 



a 



1 - &, 1 



a + 1 X — 1 



x^il-xy „ fa + b,l 

2i^i , 1 a;- 

a \ a + 1 



found for instance in [23, §8.17(ii)], the generating function 
(10) can be rewritten as 



ie-2 [ Bi{ fi- ^,1 + X 



Upon expanding this we obtain the following computation- 
ally more accessible generating function for Ls^''^^^j^ (tt): 



Theorem 1. For 2|/i| < A < 1 me have 

n,k>0 



! k\ 



+ n 



(11) 



We now show how the log-sine integrals Lsl*^' (tt) can quite 
comfortably be extracted from (11) by differentiating its 
right-hand side. The case n = is covered by: 



Proposition 1. We have 



d//'= dA" 



A = 
^ = 



-{in) 



B(m + l,k + l) 



Proof. This may be deduced from 



x-y 



- E 



m,fc>0 

E 

m,fe>0 



(fc + m + 1)! 



E B(m+l,fc + l) 



m\ k\ 



upon setting x = i-n\/2 and y = ivr^. □ 

The next proposition is most helpful in differentiation of 
the right-hand side of (11) for n > 1, Here, we denote a 
multiple harmonic number by 



Hi' 



E 



If Q = we set H, 



[0] 



n>zi >i2> ■•■>ia 

■- 1. 



ll«2 ■ 



(12) 



Proposition 2. For n > 1 



(-1)" 



(13) 



Note that, for a > 0, 
(±1)" 



E' 



ifi'!l,=Li^.{i)„(±l) 



which shows that the evaluation of the log-sine integrals will 
involve Nielsen polylogarithms at ±1, that is polylogarithms 
of the type 'Li^^^^yb{±l). 

Using the Leibniz rule coupled with Proposition 2 to dif- 
ferentiate (11) for n > 1 and Proposition 1 in the case n = 0, 
it is possible to explicitly write Lsi*' (tt) as a finite sum of 
Nielsen polylogarithms with coefficients only being rational 
multiples of powers of tt. The process is now exemplified for 
Lsf^ (tt) and Ls<^^ (tt). 

Example 2. (Lsf > (tt)) To find Lsf (tt) we differentiate 
(11) once with respect to A and twice with respect to /i. To 
simplify computation, we exploit the fact that the result will 



be real which allows us to neglect imaginary parts: 
d/i2 dA 



In the second step we were able to drop the term correspond- 
ing to n = because its contribution — i7r''/24 is purely 
imaginary as follows a priori from Proposition 2. O 



Example 3. (Ls^^' (tt)) Similarly, setting 

Li^l,...,a„ Liai,...,a„(l) " Liaj , . . . ,a,. ( - 1) 

we obtain Lsg^' (tt) as 

- Ls« (.) = ^ E ^^i^^ - H^-.) 

^ 6(1 -(-1)") n' 
rv' r}? 

= 6Li3±i,i -6Li±i +^ Ut -~\n\{i) 

= -6Li3,M(-l) + ^C(5)-i7r\(3). 

The last form is what is automatically produced by our pro- 
gram, see Example 13, and is obtained from the previous 
expression by reducing the polylogarithms as discussed in 
Section 6. <> 



The next example hints at the rapidly growing complexity 
of these integrals, especially when compared to the evalua- 
tions given in Examples 2 and 3. 

Example 4- (Lsg^' (tt)) Proceeding as before we find 
-Ls^'> (tt) = - 24Li±i,i^i +24LiJi,i -18Li±i +12Li± 
-f 3^\(3,l)-37r\(4)-h^ 



24Li3,i.i,i(-l) - 18Li5,i(-l) 



(14) 



In the first equality, the term 7r®/480 is the one correspond- 
ing to n = in (11) obtained from Proposition 1. The 
second form is again the automatically reduced output of 
our program. O 



Remark 1. From the form of (11) and (13) we find that 
the log-sine integrals Lsi'"' (vr) can be expressed in terms of vr 
and Nielsen polylogarithms at ±1. Using the duality results 
in [3, §6.3, and Example 2.4] the polylogarithms at —1 may 
be explicitly reexpressed as multiple polylogarithms at 1/2. 
Some examples are given in [6]. 

Particular cases of Theorem 1 have been considered in 
[15] where explicit formulae are given for hs^n^ (tt) where 

O 



3.4 Log-sine integrals at 2% 

As observed by Lewin [20, 7.9.8], log-sine integrals at 2ti 
are expressible in terms of zeta values only. If we proceed 
as in the case of evaluations at tt in (10) we find that the re- 
sulting integral now becomes expressible in terms of gamma 
functions: 



E 



27r / 

2 sin I j e'^''d6l 



= 27re'' 



(15) 



The special case /i = 0, in the light of (20) which gives 
Ls„ (27r) = 2Ls„ (vr), recovers (6). 

We may now extract log-sine integrals Lsl''' (27r) in a sim- 
ilar way as described in Section 3.1. 

Example 5. For instance, 

T (2) In \ 13 5 

We remark that this evaluation is incorrectly given in [20, 
(7.144)] as 77r^/30 underscoring an advantage of automated 
evaluations over tables (indeed, there are more misprints in 
[20] pointed out for instance in [9, 15]). O 



3.5 Log-sine-polylog integrals 

Motivated by the integrals LsLsCfe,ij defined in [14] we 
show that the considerations of Section 3.3 can be extended 
to more involved integrals including 



afclogn-fc-i ( 2sin^ ) Lid(e"')de 



Lsi''^ (7r;d) := 



On expressing Lid(e'^) as a series, rearranging, and applying 
Theorem 1, we obtain the following exponential generating 
function for Lsl'°' (tt; d): 



Corollary 1. For d> we have 

n,k>0 



k\ 



A 
2 



where 



(n - k)'^ ' 



(16) 



(17) 



We note for < e < tt that Li_i(e''') = -1/ (2sin|)^ 
Lio(e'*) = f cotf, while Lii(e'^) = - log (2 sin |) + 

i^, and Li2(e'«) = C(2) + f (| - tt) + i CI2 {0). 

Remark 2. Corresponding results for an arbitrary Dirich- 
let series La,d(a:) :— X]n>i o,nx"/n'' can be easily derived in 
the same fashion. Indeed, for 



Ls^*"' (tt; a, d) 



«'=log"-'=-i (2sin^)La,d(e''^)de 



one derives the exponential generating function (16) with 
H„a{X) replaced by 



- E 



[n-ky 



(18) 



This allows for Lsi'^' (tt; a, d) to be extracted for many num- 
ber theoretic functions. It does not however seem to cover 
any of the values of the LsLsc^^ij function defined in [14] 
that are not already covered by Corollary 1. O 

4. QUASIPERIODIC PROPERTIES 

As shown in [20, (7.1.24)], it follows from the periodicity 
of the integrand that, for integers m, 

Lsi'') (2m7r) - Ls^*^^ (2m7r - a) 

= E(-l)'"'(2m^)' (j) Lsi'!_7' (a) . (19) 

Based on this quasiperiodic property of the log-sine inte- 
grals, the results of Section 3.4 easily generalize to show 
that log-sine integrals at multiples of 2n evaluate in terms 
of zeta values. This is shown in Section 4.1. It then follows 
from (19) that log-sine integrals at general arguments can 
be reduced to log-sine integrals at arguments < cr < tt. 
This is discussed briefiy in Section 4.2. 

Example 6. In the case fc = 0, we have that 

Ls„ (2m7r) = 2mLs„ (tt) . (20) 

For fc = 1, specializing (19) to cr = 2rmi then yields 

Ls^^' (2m7r) = 2m^7rLs,i-i (tt) 

as is given in [20, (7.1.23)]. O 

4.1 Log-sine integrals at multiples of 2-k 

For odd k, specializing (19) to cr = 2m7r, we find 
2Lsi'=) (2m7r) = ^(-1)^~' (2m7r)^' ^^n-P (2"^7r) 

giving Lsn*"' (2m7r) in terms of lower order log-sine integrals. 

More generally, on setting cr = 27r in (19) and summing 
the resulting equations for increasing m in a telescoping fash- 
ion, we arrive at the following reduction. We will use the 
standard notation 



H 



(a) 



Y^k- 



for generalized harmonic sums. 
Theorem 2. For integers m > 0, 

Lsi'=) i2mn) = ^^(-1)''-^ (2^)^ (^^ H 



^-^^Lsi'-p{2n) 



Summarizing, we have thus shown that the generalized 
log-sine integrals at multiples of 2n may always be evaluated 
in terms of integrals at 2-k. In particular, Ls^*^' (2m7r) can 
always be evaluated in terms of zeta values by the methods 
of Section 3.4. 



4.2 Reduction of arguments 

A general (real) argument cr can be written uniquely as 
(7 — 2rmr ± ag where m > is an integer and < ag < n. 
It then follows from (19) and 



Ls« i-e) = {~ir^'Ls\:> (9) 



that Lsi'°' (a) equals 



Ls(^' (2m^) ± ^(±l)'--^(2m^)^' Lf^f (a,,) . (21) 

Since the evaluation of log-sine integrals at multiples of 2tt 
was explicitly treated in Section 4.1 this implies that the 
evaluation of log-sine integrals at general arguments a re- 
duces to the case of arguments < ct < tt. 

5. EVALUATIONS AT OTHER VALUES 

In this section we first discuss a method for evaluating 
the generalized log-sine integrals at arbitrary arguments in 
terms of Nielsen polylogarithms at related arguments. The 
gist of our technique originates with Fuchs ([12], [20, §7.10]). 
Related evaluations appear in [8] for Ls3 (r) to Lse (r) as well 
as in [9] for Ls„ (r) and Ls^^-* (r). 

We then specialize to evaluations at 7r/3 in Section 5.1. 
The polylogarithms arising in this case have been studied 
under the name of multiple Clausen and Glaisher values in 
[4]. In fact, the next result (22) with r = 7r/3 is a modified 
version of [4, Lemma 3.2]. We employ the notation 



ai,...,akj ' ai! • • ■ afc!(n — ai — . . . — afc)! 
for multinomial coefficients. 

Theorem 3. For < r < 2n, and nonnegative integers 
n, k such that n — k > 2, 



C(n-fc,{i}'=)-^ 



{-iry 



Li 



2+fe-j,{l}"-'=-2 



(e-) 



i^+'{-l) 

X 



(n-1)! ^ ^ 

^ ' r=0 m = 



(— TT Ls^ , ' At) ■ 

Proof. Starting with 



(22) 



Lifc.{i}'«(a) - Lifc,{i}"(l) = / ^^^^-Liilnifl d« 



and integrating by parts repeatedly, we obtain 

v^(— l)-* 

Life_j,{i}„ (q) - Lifc,{i}„ (1) 



j=0 



(_l)''-2 rMog^-2(z)Li{i}„+i(^) 



dz. 



(23) 



Letting a = e'^ and changing variables to z = e , as well 
as using 

(-log(l-.))" 

LHly^[z) - — , 



the right-hand side of (23) can be rewritten as 



(-1)' 



(A; -2)! (n + iy. 



'm"-' (-log(l-e-))' 



de. 



Since, for < 9 < 2-k and the principal branch of the loga- 
rithm, 



log(l 



e ) = log 



2 sin - 
2 



(24) 



this last integral can now be expanded in terms of general- 
ized log-sine integrals at r. 



c(fc,{ir)-E 



{-iry 



Life-j,{i}"(e'^) 

M)"-^ (-1)" 

(fc-2)! {n + iy ^ ^„ 



n + 1 
r 



r — Q in—0 

I ■|^g(fc + m-2) 



{r) 



(25) 



^ ■'^'^ri + fe — (r — m) 

Applying the MZV duality formula [3], we have 

c(fc,{ir) = c(n+2,{i}'=-^), 

and a change of variables yields the claim. □ 

We recall that the real and imaginary parts of the mul- 
tiple polylogarithms are Clausen and Glaisher functions as 
defined in (4) and (5). 

Example 7. Applying (22) with n = 4 and k = 1 and 
solving for Ls4^'' (r) yields 



Ls<'> (r) = 2C(3,l)-2Gl3,i(r) 



(3) 



1 T 



2rGl2,i(^) 
1 



(r) + i7r2LsW (r) 



^180- 

14^1 3 
16" + 6"" 



2r Gl2,i(r) 



1 2 2 

— TT T . 



For the last equality we used the trivial evaluation 

n 

hsf'> (r) = - — . (26) 



It appears that both Gl2,i (t) and Gls^i (r) are not reducible 
for r = 7r/2 or r = 27r/3. Here, reducible means express- 
ible in terms of multi zeta values and Glaisher functions of 
the same argument and lower weight. In the case r = 7r/3 
such reductions are possible. This is discussed in Example 9 
and illustrates how much less simple values at 27r/3 are than 
those at 7r/3. We remark, however, that Gl2,i (27r/3) is re- 
ducible to one-dimensional polylogarithmic terms [6]. In [1] 
explicit reductions for all weight four or less polylogarithms 
are given. O 

Remark 3. Lewin [20, 7.4.3] uses the special case k = n — 
2 of (22) to deduce a few small integer evaluations of the 
log-sine integrals Ls^"~^' ('""/S) in terms of classical Clausen 
functions. O 

In general, we can use (22) recursively to express the 
log-sine values LsJ,''' (t) in terms of multiple Clausen and 
Glaisher functions at r. 



Example 8. (22) with n = 5 and = 1 produces 
Ls<'' (r) = -6C(4, 1) + 6 Cls.i.i {r) + 6r Cla.i.i (r) 



3 T (3) / \ 3 

-Lss (r) ^. 

Applying (22) three more times to rewrite the remaining log- 
sine integrals produces an evaluation of Lsg^' (r) in terms of 



multi zeta values and Clausen functions at r. 



O 



As follows from the results of Section 5 each integral is a 
multivariable rational polynomial in tt as well as CI, Gl, 
and zeta values. These evaluations confirm those given in 
[9, Appendix A] for Lsa (-1), Ls4 (-1), and Ls6 (^). Less 
explicitely, the evaluations of Lss (-1) and Ls7 (|-) can be 
recovered from similar results in [15, 9] (which in part were 
obtained using PSLQ; we refer to Section 6 for how our 
analysis relies on PSLQ). 

The first presumed-irreducible value that occurs is 



5.1 Log-sine integrals at tt/s 

We now apply the general results obtained in Section 5 
to the evaluation of log-sine integrals at r = 7r/3. Accord- 
ingly, we encounter multiple polylogarithms at the basic 6- 
th root of unity uj := exp(i7r/3). Their real and imaginary 
parts satisfy various relations and reductions, studied in [4], 
which allow us to further treat the resulting evaluations. In 
general, these polylogarithms are more tractable than those 
at other values because TiJ — llP . 

Example 9. (Values at -1) Continuing Example 7 we have 
Using known reductions from [4] we get: 



Gl,,. 



and so arrive at 



324 



23 4 
:7r , 



19440 



17 4 
-TT . 



6480 



(27) 



(28) 



Lewin explicitly mentions (28) in the preface to [20] because 
of its "queer" nature which he compares to some of Landen's 
curious 18th century formulas. O 

Many more reduction besides (27) are known. In par- 
ticular, the one-dimensional Glaisher and Clausen functions 
reduce as follows [20]: 

Gl„ {2nx) = B„(a;)^", 



CI: 



2n+l 



2"'")(l-3"^")C(2n+l) 



(29) 



Here, Bn denotes the n-th Bernoulli polynomial. Further 
reductions can be derived for instance from the duality result 
[4, Theorem 4.4]. For low dimensions, we have built these 
reductions into our program, see Section 6. 

Example 10. (Values of Ls„ (7r/3)) The log-sine integrals 
at 7r/3 are evaluated by our program as follows: 

_Ls3r^u^.^ 

V3/ 108 
L..(I)4.C(3)40U(I 



LS5 - 



1543 
19440 ' 
15 



6 



Gl.,, (i 



oo 



En — 1 
fc = l 



3341 
1632960 



(30) 



47r ' — . I — )n" 

n— 1 \ n / 



The final evaluation is described in [4]. Extensive compu- 
tation suggests it is not expressible as a sum of products 
of one dimensional Glaisher and zeta values. Indeed, con- 
jectures are made in [4, §5] for the number of irreducibles 
at each depth. Related dimensional conjectures for polylogs 
are discussed in [26]. O 

5.2 Log-sine integrals at imaginary values 

The approach of Section 5 may be extended to evaluate 
log-sine integrals at imaginary arguments. In more usual 
terminology, these are log-sinh integrals 



Lshi*=' (a) := 



gk 1 n-i-fc 2 sinh- 
^ 2 



de 



(31) 



which are related to log-sine integrals by 

Lshl'^ (a) = Lsl'^ (ia) . 

We may derive a result along the lines of Theorem 3 by 
observing that equation (24) is replaced, when 6 = it for 
t > 0, by the simpler 



log(l - e ') = log 



2 sinh - 
2 



(32) 



This leads to: 



Theorem 4. For t > 0, and nonnegative integers n, k 
such that n — k > 2, 



C(n - k, {1}'=) - X] -J Li2+fe_j-,{i}„-fc-2 



fc, r 



{n-l)\ 



E 



LshL'=+'-' (t) . (33) 



Example 11. Let p := (1 + \/Z)/2 be the golden mean. 
Then, by applying Theorem 4 with n = 3 and k = 1, 



Lsh'^' (21ogp) = C(3)-ilog«p 



hh{p )-2Li2(/5 )logp. 

This may be further reduced, using Li2(p^^) = — log^ P 
and Li3(p~^ 
well-known 



and Li3(p-') = |C(3) - ^Tr^logp-f | log^ p, to yield the 



LshW (21ogp) = ic(3). 



The interest in this kind of evaluation stems from the fact 
that log-sinh integrals at 2 log p express values of alternating 
inverse binomial sums (the fact that log-sine integrals at 7r/3 
give inverse binomial sums is illustrated by Example 10 and 
(30)). In this case, 

Lsh« (21ogp) = iy i^^ii . 

n— 1 V n / 

More on this relation and generalizations can be found in 
each of [22, 16, 4, 2]. O 



6. REDUCING POLYLOGARITHMS 

The techniques described in Sections 3.3 and 5 for evalu- 
ating log-sine integrals in terms of multiple polylogarithms 
usually produce expressions that can be considerably re- 
duced as is illustrated in Examples 3, 4, and 9. Relations 
between polylogarithms have been the subject of many stud- 
ies [3, 2] with a special focus on (alternating) multiple zeta 
values [17, 13, 26] and, to a lesser extent, Clausen values [4]. 

There is a certain deal of choice in how to combine the 
various techniques that we present in order to evaluate log- 
sine integrals at certain values. The next example shows how 
this can be exploited to derive relations among the various 
polylogarithms involved. 

Example 12. For n = 5 and k = 2, specializing (19) to 
a = n and m = 1 yields 

Lsf ^ (27r) = 2 Lsf ^ (tt) - 47r Ls^^^ (vr) + Att^ Lss (tt) . 

By Example 5 we know that this evaluates as — 13/457r''. On 
the other hand, we may use the technique of Section 3.3 to 
reduce the log-sine integrals at tt. This leads to 

-87rLi3,i(l) + 127rLi4(l) - ^tt^ = -^tt^ 

5 45 

In simplified terms, we have derived the famous identity 

4 

C(3, 1) = jg^. Similarly, the case n — 6 and k = 2 leads to 
C(3, 1, 1) = |C(4, 1) + j^n^CiS) - C(5) which further reduces 

2 

to 2(^(5) — ^(^(3). As a final example, the case n — 7 and 
= 4 produces C(5, 1) = tI^ - K(3)^- ^ 

For the purpose of an implementation, we have built many 
reductions of multiple polylogarithms into our program. Be- 
sides some general rules, such as (29), the program contains 
a table of reductions at low weight for polylogarithms at the 
values 1 and —1, as well as Clausen and Glaisher functions 
at the values n/2, tv/2, and 2tt/3. These correspond to the 
polylogarithms that occur in the evaluation of the log-sine 
integrals at the special values n/3, tv/2, 2n/3, n which are of 
particular importance for applications as mentioned in the 
introduction. This table of reductions has been compiled 
using the integer relation finding algorithm PSLQ [2]. Its 
use is thus of heuristic nature (as opposed to the rest of the 
program which is working symbolically from the analytic 
results in this paper) and is therefore made optional. 



7. THE PROGRAM 

7.1 Basic usage 

As promised, we implemented^ the presented results for 
evaluating log-sine integrals for use in the computer algebra 
systems Mathematica and SAGE. The basic usage is very 
simple and illustrated in the next example for Mathematical . 

Example 13. Consider the log-sine integral Lsg^' (27r). The 
following self-explanatory code evaluates it in terms of poly- 
logarithms: 

LsToLi [Ls [5 ,2 ,2Pi]] 

This produces the output — 13/457r''' as in Example 5. As a 
second example, 

-LsToLi [Ls [5,0, Pi/3] ] 
results in the output 

1543/19440* Pi "5 - 6*G1 [{4 , 1} , Pi /3] 
which agrees with the evaluation in Example 10. Finally, 

LsToLi [Ls [5 , 1 , Pi] ] 
produces 

6*Li [{3 , 1 , 1} , -1] + (Pi-2*Zeta [3] )/4 
- (105*Zeta [5] ) /32 

as in Example 3. O 
Example 14 . Computing 

LsToLi [Ls[6,3,Pi/3]-2*Ls[6,l,Pi/3]] 

yields the value 20^4120 ^^'-^ thus automatically proves a 
result of Zucker [27]. A family of relations between log-sine 
integrals at 7r/3 generalizing the above has been established 
in [22]. O 

7.2 Implementation 

The conversion from log-sine integrals to polylogarithmic 
values demonstrated in Example 13 roughly proceeds as fol- 
lows: 

• First, the evaluation of Lsi''' (a) is reduced to the cases 
of < (J < TT and a — 2m,n as described in Section 4.2. 

• The cases cr = 2rmi are treated as in Section 3.4 and 
result in multiple zeta values. 

• The other cases o result in polylogarithmic values at 
e""" and are obtained using the results of Sections 3.3 
and 5. 

• Finally, especially in the physically relevant cases, var- 
ious reductions of the resulting polylogarithms are per- 
formed as outlined in Section 6. 

^The packages are freely available for download from 
http : //arminstraub . com/pub/log- sine- integrals 
^The interface in the case of SAGE is similar but may change 
slightly, especially as we hope to integrate our package into 
the core of SAGE. 



7.3 Numerical usage 

The program is also useful for numerical computations 
provided that it is coupled with efficient methods for eval- 
uating polylogarithms to high precision. It complements 
for instance the C++ library Isjk "for arbitrary-precision 
numeric evaluation of the generalized log-sine functions" de- 
scribed in [15]. 

Example 15. We evaluate 

— (¥)-I-...(t) 

Using specialized code^ such as [25], the right-hand side is 
readily evaluated to, for instance, two thousand digit pre- 
cision in about a minute. The first 1024 digits of the re- 
sult match the evaluation given in [15]. However, due to 
its implementation Isjk currently is restricted to log-sine 
functions Ls,?' {6) with fc < 9. O 
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